Bull. Astr. Soc. India (2011) 39,[i}j20 




Rotating Radio Transients 



E. F. Keane^* and M. A. McLaughlin^ 

'Max Planck Institiit fiir Radioastronomie, Aufdem Huge! 69, 53121 Bonn, Germany. 
-Department of Physics, West Virginia University, Morgantown, WV 26506, USA. 

Received 20 11 August 22; accepted 201 1 September 22 

Abstract. Over the past several years, it has become apparent that some radio pulsars 
demonstrate significant variability in their single pulse amplitude distributions. The 
Rotating Radio Transients (RRATs), pulsars discovered through their single, isolated 
pulses, are one of the more extreme manifestations of this variability. Nearly 70 of 
these objects have been found over the past several years in archival and new pulsar 
surveys. In this review, we describe these searches and their resulting discoveries. We 
then discuss radio timing algorithms and the spin-down properties of the 19 RRATs 
with phase-connected solutions. The spin-down parameters fall within the same range 
as other pulsars, with a tendency towards longer periods and higher magnetic fields. 
Next we describe follow-up observations at radio wavelengths. These show that there 
are periodic fluctuations in the pulse detection rates of some RRATs and that RRATs in 
general have similar spectra to other pulsars. X-ray detection has only been made for 
one RRAT, J1819-1458; observations have revealed absorption features and a bright 
X-ray nebula. Finally, we look to future telescopes and the progress that will be made 
with these in characterising and understanding the Galactic RRAT population. 
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1. Introduction 

Pulsars do not often conform to the well-known textbook model: a dependable 'cosmic light- 
house' from which we can detect a train of regularly-spaced pulses. Many are now known to 
exhibit very large variations in radio emission, and in some cases even appear to switch off for 
seconds, weeks, months or even years. To discover such malfunctioning lighthouses we cannot 
use techniques dependent on there being a reliable periodic signal. These sources are often much 
better (or only) detectable via their occasional bright bursts of emission, i.e. as 'radio transients' . 



'email: ekeane@mpi£r-bonn.mpg . de 



2 



E. F. Keane and M. A. McLaughlin 



The eRRATic pulsars discovered in this way are commonly dubbed Rotating Radio Transients, 
or 'RRATs'. 

In this review we first describe, in ^ the methods used in searches for pulsars and fast radio 
transients in general. We then discuss, in ^ the 2006 discovery of the first 1 1 RRAT sources and 
the initial implications of this. We proceed, in ^ with an overview of the numerous subsequent 
searches for radio transients which have been performed in the last 5 years. In ^we look at the 
large body of follow-up observations which have been made, primarily at radio wavelengths but 
spanning most wavelength ranges. We conclude, in ^ with a discussion of where these erratic 
pulsars fit within the overall pulsar population. 



2. Pulsar searches 



Pulsar observations have focused, in recent times, on observing frequencies of ~ 1 GHz, where re- 
ceivers with a few 100-MHz bandwidth are used routinely. The radiometer equation then tells us 
that such observing systems (on ~ 100-metre class single dish radio telescopes) can detect pulses 
of radio emission with peak flux densities of ~100 y/W^s mJy, where W^s is the pulse width in 
milliseconds. The well-known frequency-dependent dispersion a radio signal experiences as a 
result of its journey through the interstellar medium is fdeiay = 4.15(DM/vghz) ms, where the 
dispersion measure DM = J n^dl is the integrated electron density along the line of sight in units 
of pc cm"^ and vghz is the observing frequency in GHz. Therefo re, by channelising our band- 
width we can determine DMs and, with a model for ne(x, y, z) (e.g. Cordes & Lazioll2002 ) in the 
Galaxy, we can infer distances to sources. And so it is that we can easily detect millisecond pulses 
with fluxes >100 mJy and estimate luminosities. For Galactic sources the minimum brightness 
temperature associated with such transient radio signals is 4 x 10^^ K(5 jyD^p^)/(vGHzW'm,s), where 
S}y is the peak flux in Jy and D^pc is the distance in kpc. Radio pulsars have the highest bright- 
ness temperatures of all known radio transient sources, with most emitting singl e pulses in the 
range 10^^ - 10^^ K, with the giant pulses from the Crab pulsar as high as 10^*^ K dHankins et al 



20031) . Additionally, millisecond pulse widths, by causality, imply maximum source sizes of 



300 HVs km- Clearly the emission mechanism responsible for such signals involves non-thermal, 
coherent emission from compact objects. 



2.1 First discoveries 



Even though thought to exist since the 1930s ( Baade & Zwickv^ 1934 : Oppenheimer & Volkofl 
1939h neutron stars were seen as 'theoretical', i.e. difficult/ impossible to detec t. It was not 



until the serendipitous discovery of the first pulsar in 1967 dHewish et al.lll968l) that neutron 



stars were promoted to the forefront of astrophysical research. Although until this time radio 
astronomers had not observed at sufficiently high time resolution for their detection, this soon 
became commonplace. The first pulsars were discovered by visual inspection of the time series 
from large radio telescopes, at first with no c orrection for dispersion, but this was quickly in- 
corporated. The first ~30 pulsars (see § 3.4 of iLvne & Smith2004 and references therein for a 
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historical review ) were discov ered in this way, amongst which the most noteworthy is the Crab 
pulsar (Staelin & Reifenstein| [l968). The short period of 33 ms (although not measured in the 
discovery observations) solidified the neutron star interpretation (over the white dwarf interpret- 
ation) of pulsars. 



2.2 Periodicity searches 



It was quickly realised that Fourier domain searches would enable the detection of pulsars whose 
individual pulses were not bright enough to be detectable. For a train of pulses of duty cycle 
6 - W/P (i.e. width W recurring at a period of P) an improvement in sensitivity by a factor 
of V(l -<5)/<5, or of ~ l/^/6 for narrow pulses, can be made. Fourier searches, wide DM- 
range se arches, and t he use of extended bandwidths resulted in the discovery of many more 



pulsars ( Manchester e t al. 19961 2001 ). Faster time sampling allowed the discovery of milli- 



second pulsars (MSPs) (Back er et a 



systems were developed (iFreire et al 



1982[) and sophisti cated techniques for detecting binary 



.Ml^»:zi) and soptiisti c 
200ll:lRansomll200lh 



The standard pulsar search pipeline had become (see §6 of iLorimer & Krameiil2005l for an 
extensive overview): 



1. Observe at frequencies ~1 GHz, a trade off between the effects of scattering, dispersion and 
sky tem perature, which wo uld favour higher frequencies, and typical stee p pulsar spectral 



SKy tem perature, wnicn wo uia tavour nigner trequencies, ana typical stee p pulsar spectral 
indices (iMaron et al .l2000h . which would favour lower frequencies (see e.g. [Manchester et al, 



2001). 



2. Use as wide a bandwidth as possible to increase sensitivity, with as many frequency chan- 
nels as po ssible to reduce intra-channel smearing to enable detection of MSPs (see e.g. 
Keith et aL,20 10; Bates et al..,.201 1,) . 



3. Choose a DM range to search (perhaps based on the sky position and the expected max- 
imum DM in that direction) and for each DM in that range, dedisperse the data by shifting 
the frequency channels to compensate for frequency-dependent delays. The spacing of the 
DM trial values, or 'channels', is selected so that the error due to an incorrect DM is no 
smaller than the intrinsic dispersion smearing over a single frequency channel. 



4. Fourier transform these dedispersed time series and apply appropriate filters in the fre- 
quency domain to detect periodic signals. For each significant signal, fold the original 
dedispersed data at the candidate period (before tuning the DM and period values for an 
optimal detection). An acceleration search for binary pulsars adds either an extra layer to 
the nested search loop, over the range of acceleration values, or increases the complexity 
of the Fourier domain search at each DM trial. 
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2.3 Single pulse searches 

After the initial pulsar discoveries, and the switch to the 'standard' searching procedure, it was 
not standard to search the dedispersed time series in the time domain for bright individual pulses 
(with a few exceptions, e.g. Nicell999 ). However, there are a number of reasons why it is sensible 



to add such a 'single pulse search' to the pipeline. Firstly, we should exhaust all possible avenues 
to maximise the yield of our search, especially as it is computationally very inexpensive to add 
single pulse searches to the standard search setup. Secondly, it is not only quick, but effective, as 
there are many pulsars which are more significantly, or indeed only, detectable via single pulse 
searches. 



It is straightforward to show dMcLaughUn & Cordesll2003h that the ratio of signal-to-noise 



ratios of FFT and SP searches (the so-called 'intermittency ratio') is given by: 

_ (5/jV)sp _ A 5 peak 
ave 



(1) 



where A is a product of constants and is ^ 2, N is the number of pulse periods during the obser- 
vation, and 5 peak and 5 ave are the peak and average pulse flux densities, respectively. Depending 
on the pulse amplitude distribution and the observation specifications, a pulsar may be more ef- 
fectively (or only) detectable in either a single pulse (SP) or a Fourier domain (FF1[J) search. For 
longer observations there is a higher probability of a pulsar emitting a very bright single pulse 
(i.e. Speak » Save) and hence increasing the likelihood of detection in a single pulse search. 
However in the limit of high the periodicity search will become more effective. Many pulsar 
amplitude distributions have a 'sweet spot', i.e. a peak in r as a function of A^, between where 
these competing effects dominate (see lKeanell2010l for an examination of several families of dis- 
tributions). 

Single-pulse searches are quite straightforward. They involve searching each dedispersed 
time series for significant peaks, an exercise in matched filtering. Typically the dedispersed time 
series is searched for peaks, i.e. pulses of width fsamp, the sampling time, then 'decimated' in time 
(i.e. adjacent time samples are added) by a factor n and searched for peaks, i.e. pulses of width 
n X fsamp- Usually decimation is done in factors of two (with more accurate fits made later to 
obtain more precise pulse widths). For example, the Parkes Multi-beam Pulsar Survey (PMPS), 
with fsamp = 250 [xs, wa s recently searched for pulses with pulse-widths ranging from 2" to 2^ 



times this sampling time ( Keane et al.l2010l) . This is equivalent to convolving the time series with 



box-cars of width n x fsamp with n = 1, 2,4, 8, . . . , 512. Once significant events over a range of 
widths have been recorded for a range of DM values, a plot like that shown in Fig. [T]is examined 
for significant events at a constant non-zero DM. 

Further techniques are also often employed to remove signals that are due to terrestrial 



Although 'FFT' denotes Fast Fourier Transform, it is used here, as is common, as a synonym for all searches which 
rely on the time-averaged emission. 
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Figure 1. A typical diagnostic plot from a single pulse search, in this case produced using the sigproc 
software suite. The top line gives some summary information and the top three panels show, from left to 
right, the number of pulses as a function of signal-to-noise ratio, the number of pulses as a function of DM, 
and the signal-to-noise as a function of DM. The bottom panel shows all detected events (above some user- 
defined threshold, in this case 5 standard deviations above the mean) in the DM-time plane plotted as circles 
whose size is proportional to signal-to-noise. This particular observation is from the Parkes Multi-beam 
Pulsar Survey. It is 35-minutes in length and has been searched for significant pulses in over 300 trial DM 
values covering a range of - 2100 cm" ^^pc; DM values up to 600cm"^pc are shown here. Two pulsars 
(J 1840-0809 and J 1840-08 15) can be clearly seen as numerous significant events at constant DM values. 
The properties of PSRs J 1840-0809 and J 1840-08 15 are very similar, with periods, DMs, and 1400-MHz 
fluxes of 0.95 s, 350 pc cm"', and 2.3 mjy and 1.09 s, 233 pc cm"', and 1.4 mjy, respectively. However, the 
degrees of modulation of these pulsars are incredibly different. One can see by eye that PSR J 1840-08 15 
has a higher 'modulation index', i.e. Xj^(5i - S ^ycY /S^^^, and if it were more distant (equivalent to reducing 
our sensitivity) we would see just the very brightest events and the source may only be detectable in a single 
pulse search. The pulses at a DM of zero are due to RFI. 



interf erence, which can obscure as trophysical signals of interest, such as 'zero-DM subtrac- 
tion' (Eatou gh. Kean e & Lyne 20091) or coincidence matching whe n using multi-beam receivers 
dPeneva et a l. 2009; Kean e et al. teOlOtlBurke-Spolaor et alj|201 l\ h. The software packages usu- 
ally used (see McLaughlin ( 20 lib for a discussion) to search pulsar survey data are sigprocQ, 



SWINLEGIO 



iJl (an offshoot of sigproc but with many extra features) and prestcQ. All of these pack- 



ages are freely available. 



^http : //sourceforge . net/pro j ects/sigproc 

'https : //github . com/swinlegion 

''https : //github . com/scottransom/presto 
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3. RRATs: discovery 

With the aim of discover ing pulsars more easily detectable via their single pul ses, or indeed 



any tr ansient radio events. iMcLaughUn et alj ( 20061) searched the en tire PMPS (Ma nchester et al 

i lCordes 



2001 ) with the search algorithm described in lCordes & McLa ughUn'(^ 003i) and as outlined above. 



To date, the PMPS is the most successful pulsar surve y ever perfor med, with a yield of over 
800 pulsars /Manche ster et alj|200l ': Morris et al. 2002; Kramer et al. 2003; Hobbs et al.ll2004 



Faulkner et al. 2004; McLaugh hn et al.i20 06; Lorimer et al. 2006; Keith et al. 200^|Eatoughetal 



2009; Eatough 2009; Keane et al.l bOlOt fEatough et al. 2010; Keane et al. 201 1). A decade later 
this number is still growing as the survey is currently being reprocessed with a number of new 
algorithms. 



3.1 Eleven new sources 



McLaughlin et al. I (l2006h identified 1 1 new sources which were not detectable in PMPS period- 



icity searches, resounding proof of the value of single pulse searches. They were characterised 
by repeated bursts at a constant DM and were at first referred to as Repeating Radio Transients. 
However, before the initial publication, with follow-up observations of between 8 and 30 hours, 
and the detection of between 4 and 229 pulses from each source, it became clear that there were 
underlying periodicities which ranged from 1-7 seconds, typical neutron star rotation periods. 
Thence these sources were recognised to be rotating neutron stars and re-named Rotating Radio 
Transients (RRATs). Further evidence of the neutron star nature of these 1 1 sources included 
the facts that the pulse widths of a few milliseconds implied (by causality) very compact sources 
(and were close to the neutron star dynamical timescale), their inferred minimum brightness tem- 
peratures were similar to those of the known pulsars, positive period derivatives with typical 
pulsar-like values were measured for three sources. Additionally, shortly after th e discovery, one 



of the sources was seen in thermal X-rays as expected for a cooling neutron star (IRevnolds et al 
2006h . 



3.2 First ideas, questions, and implications 

The discovery of the first 1 1 RRATs in the PMPS was followed by a number of papers presenting 
theories to explain the irregularity of the detected radio emission. Radio pulses were detected at 
rates as low as once per three hours to as frequently as one every few minutes. The question of 
interest to many was whether the sources were truly 'off' (or to use pulsar parlance, 'nulling') 
during the rotation periods where no pulses were detected or whether they had some weak under- 
lying emission. The non-detection in periodicity searches of PSR J1819-1458, the source with 
the br ightest and the highest n umber of detected pulses, implied a lower limit on Speak/S^ve of 



ig nest n i 

-200 (iMcLaughlin et al]l2006l) . For J1317-5759 this hmit was a fac tor of -125. Such numb ers 



are not unheard of, e.g. in PSR B0656+14 this ratio is -420, leading l^ltevrede et aP (|2006') to 



point out that any pulsar with similar pulse-to-pulse modulation would be detected as a RRAT if 
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at a sufficient distance from Earth. If B0656+14were moved ~ 12 times further away to a distance 
of 3.5 kpc it would only be detectable in a single pulse search (and we note that the distances of 
the 1 1 PMPS RRATs ranged from ~2-l kpc). For the other nine sources there simply were not 
enough data to determine limits on S'peak/'S'ave- Nonetheless, many authors proposed different 
mechanisms for how switching 'on' and 'off' pulsar-like emission could occur, applying these 
to the 1 1 PMPS RRATs assuming they were 'off' when not detected. These ideas involved, to 
name but a few, re-activating ceased pulsar emission due to transient disturbances in supernova 
fall-back discs (Li 2 006) . the release of plasma trapped in radiation belts (Luo & Melrose 20^ 
or the interaction with proposed asteroidal material around the star (ICordes & Shannonll2008l) . 

Another implication from the discovery of the 1 1 PMPS RRATs was that there remained a 
large number of undiscovered neutron stars of this type, i.e. those with sporadic and/or highly- 
modulated radio emission which are more easily or only detectable in single-pulse searches. Con- 
sidering the likelihood of bursts being detected in the 35-minute pointings of the PMPS, the num- 
ber and distribution (e.g. in 'radio luminosit y' SD^), the extent of contamination of the survey 



due to interference and beaming effects, McLaughUn et al.l (l2006i) estimated the total number of 



such sources in the Galaxy to be larger than the 'normal' pulsars (i.e. those much more easily 
detectable in periodicity searches) by as much as a factor of four. This of course was an estimate 
based on a small number of sources and on a number of parameters whose values were not well- 
determined. However, even if over-estimated by an order of magnitude this expected number is 
still a very large fraction of the overall pulsar population. Another suggestion made in the initial 
discovery paper was that these sources represented a distinct and previously unknown population 
of bursting neutron stars. 



An estimate of the numbers and birthrates of Galactic neutron stars was made bv lKeane & Kramer 



(120081) to examine the implications of this. The conclusions were that if the neutron star classes, 
not just the 'normal' pulsai-s and the new PMPS RRATs, but also the X INS (X-ray Isolated Neut- 
ron Stars, see e.g. Kaplan 2008,]^ magnetars (Thomps on & DuncarJ 1996) . and CCOs or Central 
Compact Objects (Halper n & GotthelfibOlOl) were distinct, then the observed Galactic supernova 
rate (which is presumably greater than or equal to the neutron star birthrate) seemed to be in- 
sufficient. Identifying the populations with each other could remove this problem, leading to 
the suggestion that the PMPS RRATs and 'normal' pulsars were part of the same underlying 
population. 



The discovery of the first 1 1 PMPS RRATs raised questions regarding the pulsar emission 
mechanism(s) as well as our knowledge of how many neutron stars there are in our Galaxy. To 
address these questions a large observational campaign of the 1 1 known sources was necessary, 
as well as systematic searches to identify the large predicted number of sources. Below, in ^ we 
describe the numerous successful searches for new RRATs in various radio surveys, and then, in 
^ we describe the follow-up observation campaigns both at radio and other wavelengths. 



^These sources have also been refeired to as XDINS (X-ray Dim Isolated Neutron Stai's), INS (Isolated Neutron Stars), 
the 'Magnificent Seven', and RQNS (radio quiet neutron stars) in various publications. 
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4. The search continues 



There are several ways in which one might define a RRAT. One could argue that they are simply 
one extreme of the normal pulsar intermittency spectrum and therefore require no definition. 
However, for the purposes of this review and other papers, one could simply call RRATs those 
pulsars discovered only through a single-pulse search, as in McLaughlin et al. (2006). Or, the 
distinguishing charact eristic of a RRAT could be a pulsar from which we see predominantly 
isolated pulses (see e.g. Burke-Spolaor & Bailesl2O10l) . These are both viable but for the purposes 
of this review, we propose: 

A RRAT is a repeating radio source, with underlying periodicity, which is more significantly 
detectable via its single pulses than in periodicity searches. 

A pulsar is then labeled a 'RRAT' if, for the observing specifications concerned, r > 1 (see 
Equation [T]i- A pulsar does not have to exhibit long timescale 'nulling' (i.e. sustained peri- 
ods where radio emission is 'off') to satisfy this criterion. Clearly the definition is rather arbit- 
rary as a pulsar may be a 'RRAT' in one survey but not in another — hence the use of 'PMPS 
RRAT' etc. At the time of wr iting, acc ording to this definition, 61 pulsars have been identified 
as RRATs (" McLaughlin et al. 2006; Hessels et al. 2008; Shitov et al. 2009; Keane et al."2010!; 
Burke-Spolaor & Baile&.2010: Rubio-Herrera2010:,Keane et al..201 k.Bur ke-Spolao r et al..201 IblFI 



in various surveys, which we now review in turn. This count includes all repeating sources dis- 
covered in the above papers through single-pulse searches. 



4.1 Archival Parkes surveys 



The PMPS, wherein the original 1 1 RRAT sources were discovered, surveyed the Galactic plane 
in the longitude range 260° to 50° and latitude range \b\ < 5° at a central observing frequency 
of 1374 MHz spanned by 96 x 3 MHz channels sampled every 250 ^is. Two further surveys 
were performed with the Parkes telescope at inter mediate and high Galactic latitudes of 5° < 



< 30° (.Edwards et al.1 120011: iJacobv et al.ll2009i) in the same longitude range as the PMPS. 



The surveys used almost identical survey specific ations to the PMPS, excep t for a faster sampling 
time of 125 (xs and shorter 4.4-minute pointings. ^ Burke-Spolaor & Bailes (1201 oh analysed both 
of these surveys, performing single pulse searches and discovering 14 new sources. 



^Different numbers of RRATs are reported in several publications, for a number of reasons. Firstly, sources do not 
have a constant value of r between observations with different specifics, or even between identical observing setups 
as the sampling of the pulsars' amplitude distribution functions will always yield different 5peak/'^avc values. This is 
a natural consequence of the loose 'definition' of RRATs and we re-emphasise that the label should not be taken too 
strictly. Secondly, there have been different, even looser, 'definitions' used by diff'erent authors (the one given here 
follows [keane et aPbOl ih . none of which have been strictly applied, e.g. if sources must be completely undetectable in 
FFT searches then there will be fewer RRATs and if it is simply the discovery method which decides the label then there 
are many more pulsars which should be labeled so, e.g. the Crab pulsar. 
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With an estimated ~50% o f the expected RRAT sources in the PMPS obscured by radio 
frequency interference ([McLaughlin et al. 2006) and the development of some algorithms for 
its effective mitigation (see e.g. lEatough et al.ll2009h . iKeane et alj ( |2010[ 1201 lb re-analysed the 
P MPS, discovering 17 n ew sources. Another notable search of the Parkes archives was performed 
by iLorimer et al.l (l2007h and resulted in the discovery of the so-called 'Lorimer Burst', a 5-ms 
duration burst from an unknown apparently extra-galactic source. With many such signals expec- 
ted to be found in pulsar surveys, the DM range of single pulse searches now commonly reaches 
well beyond the 'edge of the Galaxy', i.e. the highest DM expected for a Galactic sou rce based 
on so me model for the electron density distributi on (e.g . the N E2001 model of ICordes & Lazio 
20021) . Indeed one of the bursts reported in iKeane et al.l (1201 11) is also possibly extra-galactic in 
origin. 



All of these high-impact scientific discoveries show the value of searching archival datasets 
(even those which have been examined many times already). At present there remain many 
archival Parkes surveys which have not been subjected to single pulse searches. These data, 
which are publicl y available, are a s ource with huge potential for discovery. The recent and 
ongoing efforts bv lHobbs et al.l (1201 ih to make the Parkes pulsar archives easily accessible to all 
online will no doubt help in the endeavour to exhaustively search aU available data for maximum 
scientific output. 



4.2 The HTRU 



The latest and ongoing pulsar search occurring at the Parkes Telescope is the High Time Res- 
olution Universe (HTRU, aka HTRU South) survey dKeithet al.ll2010l) . This survey is being 
performed at a central frequency of 1352 MHz with a usable (total) bandwidth of 340 (400) 
MHz. The frequency and time resolution are 390 kHz and 64 [.is respectively, ~8 and ~4 times 
higher than that of the PMPS. This has allowed the d iscovery of intrins ically narrower pulses, 
enabling more distant MSPs to be identified (see e.g. Bates et al.ll201 11) . The greater sensitiv- 
ity, coupled with the fact that the HTRU survey covers areas outs ide of the PMPS sky, means 
the expected yield of new pulsars is high, numbering around 400 ( Keith et al. 2010l) . Recently 
Burke-Spolaor et al.l ( 1201 lb() have performed a single pulse search on the HTRU data. The search 
covered 0%, 23.5% and 0.4% respectively of the HTRU low-, intermediate- and high-latitude 
pointings and resulted in the discovery of 1 1 new sources. The twin northern hemisphere HTRU 
survey will be performed with the Effelsberg Telescope in Germany and has recently commenced. 



4.3 Low-frequency Green Bank Telescope searches 



Single-pulses searches have been incorporated in several low-f requency pulsar sur veys with the 
Green Bank Telescope (GBT). The first low-frequency survey dHessels et al. 2008 ) used a cent- 
ral frequency of 350 MHz with 50-MHz bandwidth, 1024 frequency channels, and 81.92-[,is 
sampling. This low frequency, combined with the sensitivity of the GBT, makes this survey very 
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sensitive to faint, nearby pulsars. The survey covered roughly 1000 square degrees at \b\ < 5.5° 
and 75° < / < 165° and resulted in the discovery of 33 new pulsars, five of which were found 
through the single-pulse search. Two of these pulsars are very spo radic. One is undetectable 



for nearly 97% of the time and another, which 'Hessel s et al.l (120081) classify as the only RRAT 



detected in their survey, shows individual bright bursts separated by 100s - 1000s of seconds. 

A second, similar GBT survey was undertaken in the (Northern) summer and autumn of 2007. 
During this time, 10300 deg^ between declinations of -21° and +26° were covered in 'drift- 
scan' mode during the GBT track repair, with an effective integration time of 150 s. A central 
frequency of 350 MHz, 50-MHz bandwidth, 2048 frequency channels, and 81.92-|^is sampling 
time were used. A total of 26 pulsars were discovered, including one object through the single- 
pulse search. This pulsar is, however, typically detectable through an FFT. A portion of the data 
(20 Terabytes, or 1500 deg^) has bee n set aside for anal ysis by high-school students in the Pulsar 
Seai-ch Collaboratory (PSC) project dRosen et al.ll2010l) . 



Encouraged by these successes, a new 350-MHz GBT survey of the Northern celestial cap 
is currently being undertaken, with the primary goal of searching for pulsars to include in pulsar 
timing arrays for gravitational wave detection. This survey will cover declinations greater than 
38° at 350 MHz with two-minute integrations using 100-MHz bandwidth, 4096 frequency chan- 
nels, and 81.92-[xs sampling time. Given the number of sources detected in single-pulse searches 
in previous GBT surveys, this new survey is expected to discover dozens of RRATs. 



4.4 The PALFA survey 



The R\LFA (Pulsar Arecibo L-band Feed Array) survey uses a seven-beam receiver on the 305- 
m Arecibo telescope to cover the region \b\ < 5° and 40° < Z < 75°, 170° < I < 210° with 
300-MHz of bandw idth, 256 frequency channels, 64-\is sampling time, and 265-s integrations 
(ICordes et al.|l2006l) . This survey has thu s far resulted in 56 pulsar discoveries, with seven of 



those made through a single-pulse search (IDeneva et al. I l2009b . Several of these objects are very 
intermittent, with less than 10 pulses detected in many hours of follow-up observations. This 
survey should ultimately discover several hundreds of new pulsars, with tens of these identified 
as RRATs. 



4.5 Other 



In addition to the above so urces, the discovery of one RRAT source at the Puschino Observatory 



in Russia was reported b v IShitov et al. I (l2009h . Furthermore the Westerbork Synthesis Radio 



Telescope's 8gr8 survey ( Janssen et alJ 2009! ) of the Cygnus region at 328 MHz has recently 
reported the detection of four RRAT sources (lRubio-Herrera.20 1 0.) . 
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5. Follow-up work 

Many of the sources found in the surveys described above in ^have been subject to extensive 
follow-up observations, which we now describe. 



5.1 Radio 



When new pulsars are discovered, one typically does a 'gridding' observation to refine the po- 
sition (as the initial survey position can be quite uncertain, e.g. +7 arcmin for the PMPS). This 
facilitates observations at higher radio frequencies and also obtaining a timing solution. How- 
ever, the interpretation of gridding observations is difficult for RRATs (i.e. in the case of a non- 
detection it is unclear whether one was off source or whether the RRAT was 'off'). Therefore, 
refining positions must typically wait for a phase-connected timing solution, requiring a year's 
observation to remove the effects of the Earth's motion around the Sun. While most pulsars are 
timed through their average profiles, RRATs are timed from their isolated pulses and often a 
large amount of time may be required to accumulate enough pulses to achieve a solution. Con- 
sequently, some RRATs (especially those with the lowest detected burst rates) still do not have 
solutions after several years and some are so sporadic that it is not possible to determine a timing 
solution. When solutions are determined, errors on timing parameters can be much higher than 
for nor mal pulsars due to pulse-to-pulse jitter and varying numbers of components in the single 
pulses (iLvne et alJl2009HHu et al.ll201 ih . 



Through observations with the Parkes, Lovell, Green Bank , and Arecibo telescopes, tim- 
ing solutions have been determined now for 19 PMP S RRATs dMcLaughlin et all 120061 120091: 
Lvne et alj2009l: Keane et al]|201 1 ; Bovles et al.ll201 1 ). Achieving phase-connected solutions al- 
lows us to place these sources on the period-period derivative diagram, a commonly used scatter 
plot for pulsar classifications, as in Fig. [S] These timing solutions have shown that there is no 
evidence of any binary companions for these sources. As can be seen from Fig. |2] while the 
spin-down properties of RRATs and other pulsars overlap, RRATs in general have longer periods 
and higher magnetic fields than other pulsars. A Kolmogorov-Smirnov test shows probabilities 
of 1 X 10 and 2 X 10"^ that the periods and magnetic fields of RRATs are drawn from the 
same popula tion as normal non - recycl ed pulsars, respectively. These probabilities are lower than 
calculated in McLaughlin et al. ( 20091) from a sample of nine RRATs, indicating that this trend is 
robust. However, it is unclear how much of this is intrinsic and how much is due to the selection 
effects of surveys. One interesting observation is that, while many RRATs populate the same re- 
gion of P-P space as other pulsars, several RRATs occupy P-P space devoid of normal pulsars 
and closer to the population of XINS. This cou ld imply an evolutionary relationship between 
these objects, as suggested bv lPopov et al. I (l2006h . 



Timing solutions have also facilitated high-energy follow-up observations (discussed below), 
and enabled monitoring for timin g irregularities. I n particular, they allowed the detection of two 
glitches from PSR J1819-1458 dLvneet al.ll2009l) . The glitch recovery is anomalous in that it 
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Figure 2. A portion of the pulsar P — P diagram showing all known pulsars. The 19 sources (18 PMPS 
'RRATs' lMcLaughlinetal.ll2009l: iKeane et alJboul and one GBT 350-MHz RRAT iBovles et al .11201 ih 
which have determined period derivatives are shown as open black stars. Magnetars are marked as plus 
signs and XINS as crosses. Also over-plotted are lines of constant 'characteristic age' and vacuum dipole 
estimate of the magnetic field (see e.g. lLorimer & KrameJ |200^ for definitions of all of these quantities). 
The periods of the RRATs with no measured period derivatives are shown at the top of the diagram. 



results in a net decrease in the st ar's slow-down ra te, the opposite trend to that in all other pulsar 
glitches (lEspinoza et alJl201 Ih . iLvne et al. I (l2009h speculate that if this behaviour were typical 
for PSR J1819-1458 it may have evolved from the magnetar region of the P - P diagram. This 
source is being moni tored regularly with the Loyell te lescope in the UK and with the Urumqi 
telescope in China ( Esamdin et al ] |2008l:lHuet"al]|201li) . No further ghtches have yet been seen. 
No glitches have been detected as yet from any of the other so urces. PSR J1819-1458 is also the 
only source which has measured radio polarisation properties. iKarastergiou et al. I (l2009h found it 
to be similar to other slow, low spindown-energy pulsars and is seen to have a 'normal' S-shaped 
swing of the polarisation position angle. 



A recent analysis by iPalhvaguru et al. 



mm of eight of the original 1 1 PMPS RRATs 
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searched for short (approximately hours) and long-term (weeks to years) periodicities in the 
detected pulse arrival times. They found significant periodicities in the pulse arrival times for 
six RRATs with timescales ranging from minutes to years and were not found in sequences of 
simulated random pulses. They also found that for the two brightest and most prolific sources 
(J1317-5759 and PSR J1819-1458) there were more instances of consecutive pulses seen than 
would be expected from a random distribution of pulse times (though most pulses, even for these 
sources, are s olitary). An o ther a nalysis of RRAT spectra based on multi-frequency radio ob- 
servations by iMiller et al. ( 201 1 ) showed that RRATs have typical steep pulsar spectra with 
spe ctral indices ~1.7 (iMiller et al. 201 lb unlike the flat spectra seen in radio-emitting magnet- 
ars (ICamilo et al ] l2006l] l2008h . \\^ note, however, that these spectral indices are based on meas- 
urements at frequencies >700 MHz and that lower-frequency spectral turnovers are possible. 



5.2 Optical and Infrared 



PSR J1819-1458 has been obse rved in the optical in the search o f a counterpart to the radio 
and X-ray (see below) source by Dhillo n. Marsh & Litt lefair ('2006*). They observed this source 
with ULTRACAM (Dhillon et al . 2007) on the 4.2-metre William Herschel Telescope, acquiring 
97100 images, each with an ~18 ms exposure time. No detection of any optical bursts were 
made with a limiting burst magnitude of /' - 16.6. To improve upon this, a follow-up study 
again using the William Herschel Telescope as well as on the New Technology Telescope was 
performed in tande m with simultaneous radio observations with the Lovell Telescope at Jodrell 



Bank (Dhillon e t al.ll201 Ih . The goal was to identify the times when radio pulses were detected 
and examine the corresponding optical frames. This is the best way to obtain an optical limit un- 
der the assumption t hat the optical emis sion is also transient and connected to the radio emission. 



as seen in the Crab (' Shearer et al.ll2003 ). No evidence for optical bursts or increased optical flux 



during radio bursts was found. Besides these observations of PSR J1819-1458 there have been 
no other reported optical observations of RRATs. 



Rea et al. I (l2010h attempted to detect infrared counterparts of the first three RRATs with tim- 



ing solutions, and hence the first with accurately determined positions. Using the Very Large 
Telescope in Chile they found no infrared counterpart to J13 17-5759 down to a limiting mag- 
nitude of Ks ~ 21. For J1913+1330 the observations used an incorrect source position due to a 
slight (yet crucial) inaccuracy in the originally published timing solution, so that the source lay 
outside of their field of view. A putative counterpart is reported for PSR J1819-1458 on the edge 
of the 1-cr error circle of the X-ray position, but .Rea et al.. (.2010.) conclude that the probability of 
a chance alignment is high. 



5.3 X-ray 



Se veral X-ray obse r vation s of PSR J18 19- 1458 have been made, at first a serendipitous detection 
by [Reynolds et al.l (120061) using Chandra (30 ks) and then a targeted observation with XMM 
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Newton (43 ks) bv lMcLaughlin et al. These revealed a thermal spectrum with kT~ 14-0 eV, 

characteristic of a cooling neutron star Pulsations at the 4.2-s radio-determined period were 
also observed, as well as a 0.5 keV spectral absorption feature, tentatively identified with proton 
cyclotron resonant scattering. Such features enable an independent estimate of the magnetic field. 
The inferred magnetic field strength from the cyclotron interpretation is 2 x 10'"* G. Equating this 
with the (commonly used but incorrect) vacuum dipole estimate for the magnetic field predicts the 
magnetic and rotational axes to be misaligned b y a - 15°. Interest ingly, equating the cyclotron 
estimate with that from the force-free model of ISpitkovskv (l2006t) yields sin^ a < 0, meaning 
either the cyclotron interpretation or the force-free magnetic field strength esti mate is incorrect . 
The spectral feature is confirmed in Chandra observations (30 ks) made by Rea et al. ' ('20091) 
which also reveal evidence for extended emission up to 5.5 arcsec from the star. Rea et al. (2009|) 
suggest that this emission is due to a pulsar wind nebula but c annot explain the power source 
given what would be a large required X-ray efficiency (see e.g. iGaensler & S land 120061) . They 
suggest that the PWN could be powered by the large magnetic field of this RRAT. A detailed 
study of PSR J1819-1458 using XMM Newton (94 ks) in tandem \yith th ree radio telescopes has 
reported initial results confirming the earlier findings d MiUer et al.ll201 11) . with a comprehensive 
analysis, including a correlation analysis of radio/X-ray flux, in preparation. 



In addition to PS R J1819-1458, PSR J1317-5759 was observe d with XMM-Newton (32 ks) 
but with no detection dRea & McLaughlinll2008l:lMcLaughlinll2009l) . PSR J1 913-H1330 was also 
observed, using Swift (9.3 ks), but again w ith an inaccurate source position ( Rea & McLaughhnI 



20081: iMcLaughlinI 120091: iRea et al.ll2010(). When im proved astrometry became available for 



PSR J0847-43 1 6 and PSR J 1 846-0257 (IMcLaughUn et al. 2009), these were both observed with 
Chandra (10.7 ks and 196.1 ks r espectively) but neith er were detected, with the corresponding 
limits obtained being reported in 'Kaplan et al.' ^20091). An X-ray candidate for PSR J1911+00 
has also been suggested by Hoffman et al. (2006) which fortuitously occupies the same field as 
Aql X-1, a much-observed low-mass X-ray binary, but no pulsations were detected and the spec- 
trum is not thermal. With by far the lowest rate of detected bursts of the original PMPS RRATs, 
PSR J 19 11 +00 is unlikely to ever have an accurate position determined, and hence progress in 
identifying this X-ray source is unlikely. 



5.4 Gamma-ray 



With the advent of the Fermi satellite, we now know of nearly 100 rotation-powered pulsars that 
are detectable at gamma-ray energies. The gamma-ray luminosity is roughly proportional to the 
spin-down luminosity, with all of the Fermi sources in the first catalog h aving Es grea ter than 



10^"* erg s ' . Only one RRAT, J 1554-5209, has E greater than 10^"^ erg s ' ; 



Keanel(l20 id) folded 



Fermi photons using the radio ephemeris and found no evidence for gamma-ray emission from 
this source. The remaining 1 8 RRATs hav e Es less than 10^^ erg s ' and are not predicted to be 
gamma-ray sources, though Keand ( 2010l) detect a folded profile of 4-cr significance with Fermi 
for J1819-1458. However, they only detect emission in half of the Fermi dataset. A longer 
timespan of Fermi observations is necessary to gauge the reality of this signal. 
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6. Where do these sources fit in? 

With the discovery of so many transient radio neutron stars in single pulse searched we now 
know of RRATs covering a large range in g, the fraction of the observed pulse periods where we 
do not detect a pulse (so that 1 - g' is an upper limit on the nulling fract ion). The p ulsars identified 
as RRATs span the range 10^^ ^ g ^ 10 ' and, as noted by Burke- Spolaor et a l. (2011b), if we 
assume these sources are in deed nulling, then they seem to be an extension of the nulling pulsars 



discovered in the PMPS bv lWang. Manchester & Johnston! (l2007h 



As well as in nulling fraction, the distribution seems to be continuous in the typical 'off' 
timesc ales for these 'nullers' and 'RRATs', spanning the range 10 s < foif < 10"* s dBurke-Spolaor et al 



201 Ibl) . This hints at a continuum in intermittency timescales, but beyond hour-long timescales 



there is a dearth of known objects until we reach the so-called 'intermittent pulsars' , PSRs B 1931-H24, 



J 1 832-1-0029, and J 1 107-5907 where the relevant timescales range from ~ 1 week to 1 year (IKramer et al, 



2006l:lKramei]l2008h . Of course the lack of objects with such long 'off' timescales is to be expec 



ted due to the difficulty in discovering (and confirmingfl) such sources. 



6.1 The trouble with timescales 



The timescales on which pulsars exhibit variability span more than 16 orders of magnitude, from 
the nano-structure seen in the Crab pulsar (~10"'^ s) to the approximately year-long timescales 
of intermittent pulsars (~10^ s). However it is very difficult to explain such a large range of 
timescales. It is an open question as to how stable pulsar emission can occur (at all!) and for days 
or weeks before apparently switching off. But the picture is even more complicated as periodic 
switching back and forth between (at least two) stable magnetospheric states has been reported 
in six pulsars studied by ' Lvne et aL I (l2010h . They report correlated period derivative and pulse 



shape changes which, when understood to be the result of a switch between two stable states, can 
remove the so-called 'timing noise' in these pulsars. Such results hint that observed phenomena 
like nulling (or the more general phenomena of mode-changing, and therefore RRATs), pulse- 
shape changes and timing noise are all one and the same and a result of magnetospheres switching 
between stable states. 

A number of stable solutions have been obtained n umerically for force-free magnetospheres 
dContopoulos. Kazanas & Fendtlll999l:lTimokhinll2006l) . Switching between these stable config 



urations is possible if there is a sudden depletion of charges in the mag netosphere. IContopoulos 



(l2005h iTimokhin (20l3), 

and lRosen. McLaughhn & Thompsonl(l201lh have even discussed how 



these switches would alter the radio emission and result in observed moding/nulling, but no 



'For a complete list of the pulsars which have been found in single pulse searches since 2006 refer to the website 
http : //www . as . wvu . edu/~pulsar/rr atalog/. 

'*In general a pulsai' is not considered to be a truly new discovery until it has been confirmed in at least one follow-up 
observation, post-survey. 
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trigger mechanism for such an event is known. Perhaps the most puzzling aspect of all is that 
the mechanism for these switches is quasi-periodic in nature. One possible cause of quasi- 
periodicities cou ld be due to modulations o f magnetospheric particle density by asteroids orbiting 
the neutron star ( ICordes & Shannonll2008 l). However, due to the low mass needed in these belts, 
they may be undetectable through reflected radio emission until the SKA or directly without a 
larger infrared telescope. 



6.2 The population of radio transients 



We note that the large projected population of pulsars likely to be detected as RRATs may not 
be a cause for concern. This is because it seems that many can be accounted for as weak/distant 
sources with high modulation indices. If the undetected underlying emission from such pulsars 
is above the typical minimum radio luminosity above which the pulsar population is estimated 
(Lmin ~ 0.1 mJykpc^) then these pulsars are akeady accounted for within the scaling factors used 
to calculate the estimated total number of pul sars (above Lmin) from the observed sample (see e.g. 
Lorimer et a l.'20Q6V Ridlev & Lorimej20 10'). However, if some of the RRATs are indeed nulling 



pulsars, or pulsars whose underlying emissi on is below the arbitrarily-imposed luminosity cut 



off' (see e.g. lFaucher-Giguere & Kaspi l2006h . then they will indeed contribute to the number of 



Galactic pulsars. Pulsar population estimates do not take account of nulling when determining 
the likelihood of detection of pulsars of various types (which consequently inform how known 
selection effects are removed to reach an estimate for the total population), another contribution 
(in addition to a luminosity-limited calculation) to the fact that the total number of active pulsars 
in the Galaxy should be larger than estimated through these methods. 

The search for radio transients has also resulted in the discovery of other interesting sources. 
Mos t notably, as we have mentioned abov e, is the Lorimer Burst, of which there may now be 



two (iLorimer et al .bOOTtlKeane et alJ201 l|). As radio frequency interference is a major hindrance 
to such searches, it is not surprising that previously unknown types of terrestrial signals have also 
been identified in transient searches. Burke-Spolaor et al. (201 la) report the detection of signals, 
referred to as 'perytons', with anomalous dispersion characteristics which are not of astrophysical 
origir@. Such terrestrial signals are a good illustration of the limits of single-dish radio telescope 
surveys. This is another argument in favour of the move towards radio arrays where signals from 
many elements (not necessarily dishes) can be correlated as appropriate to identify the source of 
terrestrial signals and to remove them. 



They are detected in all 13 beams of the Parkes 21 -cm multi-beam receiver, which, as iBurke-Spolaor et alj J20TT^ 
note, is impossible for a boresight signal. 
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6.3 RRATs with TNG instruments 



Many of the next generation of radio telescopes which will be used for pulsar and fast tran- 
sient searches are array instruments. The potential of transient searches has now truly been 
recognised, with all of these major telescopes (e.g. LOFAR, FAST, ASKAP, Mee r KAT, SKA) 



having pulsars and fast transients as key science objectives dJohnston et al.ll2008t iBooth et al 



2009 : van LeeuwenI 201 1 : Nan et al. 2011). LOFAR, for instance, curren tly in its comrnission - 



ing stage, but is already providing scientific results in the field of puls ars ( Stapp ers et al.''201 
LOF AR predicts the discovery of ~1000 pulsars in an all-sky survey (Ivan Lee uwen & Stapper 



l2010i) : possibly hundreds of RRATs will be discovered. Even greater advances will come with 
surveys for transients in interferometric images on short time scales. They will enable accur- 
ate positional localisation at the time of discovery, facilitating timing solutions and immediate 
high-energy follow-up. As a result of the transition to 'software telescopes' like LOFAR, much 
software development is needed. Already, much work is underw ay in, e.g . live t ransient searches 
using graphical processor units as recently demonstrated by Mag ro et al. ( 201 ih . Substantial ad- 
vancements like this are essential to enable us to deal with data from the SKA by the end of the 
decade. 
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